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Abstract 

We search for rare FCNC charm decays of the form — » h + l + l'~ , where A+ is a charm hadron, 
h is a pion, kaon or proton, and is an electron or a muon. In the pion and kaon modes, we 
study both D + and Df decays, while in the proton modes we study decays. Based on a data 
sample of 288 fb" 1 of e + e~ collisions collected by BABAR, we set preliminary 90% confidence level 
limits between 4 to 40xl0 -6 for the branching fractions of the different decay modes. For most 
decay modes, our analysis provides a significant improvement over previous results. 
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Figure 1: Standard model short-distance contributions to the c — > u£ + £ transition. 

1 INTRODUCTION 

In the Standard Model (SM), flavor-changing neutral-current (FCNC) processes cannot occur at 
the tree level. FCNC processes therefore provide an excellent tool for investigating the quantum 
corrections in the SM as a way to search for evidence of physics beyond the SM. FCNC processes 
have been studied extensively for K and B mesons in K° — K° and B° — B° mixing processes and in 
rare FCNC decays, such as s — > d£ + £~ , b — ► 57 and b — » s£ + £~ decays. The present measurements 
of these processes agree with SM predictions [1], but there are strong ongoing efforts to improve 
both the measurements and the theoretical predictions, and to measure new effects, such as CP 
violation, in FCNC processes. 

FCNC processes in the charm sector have received less attention and the experimental upper 
limits are currently above the SM predictions. In the SM very small signals are expected, as 
a consequence of effective Glashow-Iliopoulos-Maiani (GIM) cancellation. For instance, the c — > 
u£ + £~ transitions illustrated in Fig.^lead to branching fractions for D — > X u £ + £~ of O(10~ 8 ) [2,3]. 
This contribution is masked by the presence of long-distance contributions from intermediate vector 
resonances such as D — > X U V, V — ► £ + £~. These are predicted to have branching fractions of 
O(10 -6 ) [2,3]. In c -> u£ + £~~ transitions, the effect of these resonances can be separated from 
short-distance contributions by studying the invariant mass of the £ + £~ pair. In radiative charm 
decays, c — > wy, the long-distance contributions make it impossible to study the underlying short- 
distance physics [4]. 

Several extensions to the SM have been studied and their impact on D — > X u £ + £~ decay rates 
estimated [2,5]. The largest possible effect is expected in R-parity violating supersymmetric models. 
Depending on the size of the R-parity violating couplings, branching fractions of up to O(10 -5 ) for 
different D — > X u £ + £~ decays are possible. This is within the reach of our experimental sensitivity 
of O(10 _6 -10 -5 ), depending on the decay mode. 

There is a large group of possible FCNC charm decays to be measured. The best existing limits 
are for the branching fractions BF(D° — * £ + £'~) [6] and are O(10 -6 ) at 90% CL. In this analysis 
we search for FCNC charm decays of the form X+ — * h + £ + £'~, where the two leptons £ + and £'~ 
can each be either an electron or muon. The charge-conjugate decay modes are implied here and 
throughout this note. Current upper limits [7-10] range from no limits in some of the baryon modes 
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Table 1: Branching fractions [11] for the charm decays used for normalization. 



Decay mode Branching Fraction 

D + -> (6.2±0.6) x 10" 3 

D+ -» (36±9) x 10~ 3 

A+ -> pK-ir + (50±13) x 10~ 3 



to 5 x 1CT 6 . The FCNC decay combinations we examine are D + -> ir + £ + £'~, Df -> and 
71+ — > p£ + £'~ . Decays where the two leptons are of different flavor are lepton-family violating and 
therefore forbidden in the SM. We also search for D + — > K + £ + £'~ and — > 7r + ^ + £ /_ decays, but 
these require both quarks in the charm meson to change flavor. 

The only long-distance contributions relevant at the current experimental sensitivity are from 
D + — > 7r + and D+ — > 7r + (/> decays. The branching fractions to ir + £ + £~ through these two 
resonance decays are 1.8 x 10~ 6 and 1.1 x 10 -5 , respectively [11]. In this analysis, we measure the 
total rate of decay excluding a region around the (ft resonance in the invariant mass of the £ + £~ 
pair. 

The measured FCNC decay yields are converted into branching ratios by normalizing them to 
the yields of known charm decays. We choose normalization modes that have kinematics similar 
to the FCNC decays, so that most of the systematic effects not related to particle identification 
cancel in the branching ratio. For the D + and decays we use decays to n + (p and for the A£ 
we use A£ — > pK~ir + decays. The measured branching fractions for these modes are listed in 
Table ^ The <p decays are reconstructed in only the K + K~ decay modes. This introduces an 
additional branching fraction for <fi — > K + K~ , which is 0.491 ±0.006 [11]. We use the abbreviation 
— > tt + 4>kk to denote the Df-. — * tt + 4>, <j> — * K + K decays. 

2 THE BABAR DETECTOR AND DATASET 

The measurements are performed using data collected by the BABAR detector [12] at the PEP- 
II storage ring at SLAC. The data sample used comprises an integrated luminosity of 263 fb -1 
collected from e + e _ collisions at the T(4S) resonance and 25 fb" 1 collected 40 MeV below the T(4S) 
resonance. For event simulation we use the Monte Carlo (MC) generator EVTGEN [13] with a 
full detector simulation based on GEANT4 [14]. Signal and the A£ — > pK^n + MC events are 
generated with a 3-body phase-space distribution while Df^. — > tt + 4>kk MC events are generated 
with a Breit-Wigner for the (j) decay. All signal events are simulated as cc continuum events. 
Samples of simulated generic cc and uds continuum events and BB decays corresponding to 1.4-5 
times the recorded data sample are used to study background contributions. 

3 ANALYSIS METHOD 

Initial charm hadron candidates are formed from one track identified as either a pion, kaon or 
proton and two tracks each of which is identified as an electron or a muon. Typical electron 
(muon) identification efficiency is about 94% (60%). Hadron identification efficiencies average 
about 98%, 87% and 80% for pions, kaons and protons, respectively. The two tracks identified 
as leptons are required to have opposite charge. For candidates with a pion or kaon track, the 
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invariant mass of the charm candidate is required to be between 1.7 and 2.1 GeV/c 2 , while for the 
candidates with protons it is required to be between 2.2 and 2.4GeV/c 2 . Charm hadrons from 
continuum production typically carry most of the energy of the initial quarks. We therefore require 
the momentum, p*, of the charm hadron candidate in the e + e~ center-of-mass frame to be larger 
than 2.8 to 3.5 GeV/c depending on decay mode. This removes a large fraction of the background. 
Charm hadrons produced in B decays are kinematically limited to be below about 2.2 GeV/c and 
are therefore not used in this analysis. 

After the initial event selection, significant combinatoric background contributions remain from 
semileptonic B decays and low-multiplicity QED events. These background sources have been 
studied using events from invariant mass sidebands in data and from MC samples. The final event 
selection criteria are chosen based on these studies to minimize the expected upper limit on the 

— > h + £ + l'~ branching ratio under the assumption that no signal is present. 

The QED events are mainly radiative Bhabha, initial state radiation and two-photon events. 
These are easily identified by their low multiplicity and/or highly jet-like structure. They are 
suppressed by requiring at least five tracks in an event and a minimum event sphericity [15] of 0.13 
calculated in the e + e~ center-of-mass frame. 

We suppress the background from semileptonic B decays by rejecting events with evidence of 
neutrinos and requiring the two leptons to come from a common point. The latter is achieved by 
a tight requirement on the probability of the vertex fit x 2 {P(x 2 ) > 0.05) and on the distance of 
closest approach between the two leptons (< 250 /iin). 

The neutrinos are not directly detected and therefore show up as missing energy in the event. 
This is measured in two ways. We first calculate the total energy in an event from all reconstructed 
neutral clusters and tracks in an event, assuming the pion mass for tracks. Second we calculate the 
net transverse momentum in an event (with respect to the beam axis) by adding the momentum 
vectors of all neutral clusters and tracks. The neutral clusters are assumed to be photons. The 
more the transverse momentum deviates from zero the more likely it is to be a semileptonic B 
decay. Figure [2] shows the distributions of the two variables for D + — > 7r + e + e~ signal MC events 
and for 7r + e + e~ candidates in generic MC B + B~ events. Selecting on a linear combination of the 
two variables as indicated in the figure, 74% of the signal events are kept while 84% of the B + B~ 
events are rejected. 

In the e + e~ decay modes, there is a significant contribution at low e + e~ mass from ir° decays. 
Some of these are from decays directly to e + e~(^f) while others are tt° — » 77, where one of the 
photons converts to e + e~ . These are removed together with 7 conversions by requiring m(e + e _ ) > 
200 MeV/c 2 . 

For Df — * tt + £ + £'~ decays, we furthermore require the presence of a photon candidate with 
an energy above 100 MeV and \m(M + £'-^) - m(hl + t'~) - 0.1438 MeV/c 2 | < 0.015 MeV/c 2 . This 
selects Df mesons coming from D* + — > Dfj decays and helps reduce the non-D+ background. 

The requirement on the momentum p* of the charm hadron is decay mode dependent: 

• For D + h£+i'-, p* > 3.3GeV/c. 

• For D+ -► h£+l'-, p* > 3.5 GeV/c. 

• For A+ -> p£ + £'-, p* > 2.8GeV/c. 

For the Df^ — > it + l + £~ decay modes, we exclude events with 0.95 < m(e + e~) < 1.05GeV/c 2 
and 0.99 < m(> + //~) < l-05GeV/c 2 . The excluded regions for the two decay modes are different 
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Figure 2: Total energy vs net transverse momentum for D + — > 7r + e + e~ MC signal events (points) 
and background events from generic B + B~ events (shaded histogram) for the ir + e + e~ mode. The 
events below the red line are rejected in the final event selection. 



due to the larger radiative tails in the 7r + e + e decay mode. In order to make further cross- 
checks, we also select a sample with — > ttc/), (j) ~~ * £ + £~ candidates. For this sample we require 
0.995 GeV/c 2 < m(e + e~) < 1.030 GeV/c 2 , 1.005 GeV/c 2 < m{n + n~) < 1-030 GeV/c 2 and p* > 
3.1GeV/c. 

For the normalization decay modes the same selection is applied except for the lepton particle 
identification requirements and the restriction on the distance of closest approach between the two 
lepton candidates. Instead the kaons and pions in the — > tt + 4>kk and A£ — > pK~ir + decay 
modes are required to be specifically identified. For the Df^ — > tt + 4>kk decay modes, we further 
require the invariant mass of the kaon pair to be within 15 MeV/ c 2 of the world-average value for 
the 4> mass [11]. Since the selection of — > tt + £ + £'~ decays requires the presence of a photon from 
D* + — > Dfry decays whereas the Df — > K + i + i'~ selection does not, the selection of Df — > tt + 4>kk 
decays is performed separately for the two decay modes. 

The invariant mass distributions of the h + l + £'~ candidates are fitted using an extended un- 
binned maximum likelihood fit. The probability density function (PDF) for signal events is given 
by the so-called Crystal Ball function [16] in order to account for radiative tails: 



P CB (m;n,a,a,n) = I 2 / \" (1) 

e u -° 3 , r if m < a — aa. 



e K a > if m > jx — aa. 

na 

na— a 2 a— a(m— /t) 

The four parameters, fi, a, a and n, are obtained from fits to signal MC and kept fixed during fits 
to the data, leaving only the overall normalization as a free parameter. The fitted width of the 
Gaussian component (a) is found to lie between 6 and 9MeV/c 2 depending on the decay mode. 
The MC events have been corrected to reproduce particle identification efficiencies measured in 
various control modes. The yields from the MC fits are used to calculate the signal efficiencies, 
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which range between 0.3 and 5.3%. The modes with low efficiency are decays with two muons 
and decays requiring a photon from D* + decays. The muon identification efficiency is very low for 
muons with momentum below 1 GeV/c. 

For the normalization decay modes, the radiative effects are negligible and we use the sum of 
two Gaussian distributions with a common mean to describe the D + , and A£ signals. All 
parameters are free in the fits to data. 

The invariant mass distribution of the combinatoric background events for the signal modes are 
described by first-order polynomials. For the normalization decay modes, a second-order polynomial 
is used. The background parameters are allowed to vary freely in all cases. 

An additional background comes from hadronic charm decays where two hadrons are misiden- 
tified as leptons. Pions have a probability of about 2% (0.1%) to be identified as a muon (electron) 
depending on the pion momentum and angle. This background component is negligible in the 
signal modes with electrons and is therefore only included in decay modes with two muons. The 
shape of this background is obtained from MC samples of hadronic three-body charm decays. Each 
event is weighted according to the probability of misidentifying a pion as a muon. The misidentifi- 
cation probability is measured from data using samples of D° — > K~ir + decays. The misidentified 
hadronic charm decays are reconstructed at slightly lower h + fi + fi~ mass than the signal events. 
The peak mass is shifted by about 15 MeV/c 2 which is sufficient separation for the yield to be 
determined by the likelihood fit instead of relying on the MC prediction. 

For the charm meson modes two signals are fitted simultaneously since the D + and Df mesons 
can decay to the same final state. However, since the optimal event selection criteria are different 
for the D + and Df modes, only one of the fitted yields is used from each fit. 

For most modes we do not see a significant signal; in addition to determining a central value 
for the yields, we set upper limits on the branching ratios at 90% confidence level (CL). This is 
chosen as the point where the negative log likelihood is 1.355 above its minimum value. In order to 
always have a physical (positive) upper limit, we only consider the minimum at or above signal 
events. This generally results in conservative confidence intervals. Systematic uncertainties from 
the signal efficiency and the normalization mode, detailed in the next section, are included as a 
Gaussian constraint in the likelihood expression. 

4 SYSTEMATIC STUDIES 

Most systematic effects are expected to cancel in the branching ratio since they affect the signal 
and normalization modes equally. We therefore only have to account for differences in selection, 
acceptance and decay kinematics. Table 121 gives a summary of all the systematic uncertainties 
related to the branching ratio calculation. An additional systematic is assigned for the estimation 
of the signal yield. The details of the different uncertainties are given below. 

Systematic uncertainties related to the signal PDF parameters obtained from MC are inves- 
tigated in two ways. First, the PDF parameters for data and MC are compared in the fits to 
the normalization modes. Differences can be due either to general data-MC tracking differences 
or to uncertainty in the PDG mass used in the simulation. Second, fits to the invariant mass of 
J/ip — > £ + £~ candidates from inclusive B decays are compared between data and MC. The second 
comparison is sensitive to effects associated with lepton reconstruction. Based on these studies, 
fits with the mean mass shifted by up to 2.5 MeV/c 2 , depending on the decay mode, are performed. 
Based on the same studies, the widths of the signal PDFs are changed by ±3%. The fit giving the 
highest upper limit on the branching ratio is used for the final result. 
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Table 2: Summary of the multiplicative systematic uncertainties for all the decay modes. DOCA 
is the uncertainty due to the distance of closest approach requirement and PID is the uncertainty 
from the particle identification. 

Normalization 

Deca y mode Mode MC stat. DOCA PID p* Total 

D+ 
D+ 
D+ 
D+ 
Df 
Df 
Df 
Df 
D+ 
D+ 
D+ 
D+ 
Df 
Df 
Df 
Df 
A+- 
Af- 
At- 
Af- 
D+ 
D+ 
Df 
Df 



1 1 


3.7% 


1.0% 


0.3% 


3% 


2% 


5.3% 


1 1 


3.7% 


1.8% 


0.3% 


5% 


2% 


6.8% 


\ | 


3.7% 


1.3% 


0.3% 


3% 


2% 


5.3% 


1 1 


3.7% 


1.3% 


0.3% 


3% 


2% 


5.3% 


■ vr+e+e" 


4.0% 


1.8% 


0.3% 


3% 


2% 


5.7% 


+ + — 


A nO/ 

4.0/o 


O 00/ 

3.37o 


0.3% 


5% 


oo/ 
2% 


7.57o 


■ 7r + e + [i~ 


4.0% 


2.3% 


0.3% 


3% 


2% 


5.9% 


■ 7r + [i + e~ 


4.0% 


2.3% 


0.3% 


3% 


2% 


5.9% 


■ K+e + e- 


3.7% 


1.1% 


0.3% 


4% 


2% 


5.9% 




3.7% 


2.2% 


0.3% 


5% 


2% 


6.9% 


K+e + n~ 


3.7% 


1.5% 


0.3% 


4% 


2% 


6.0% 


K+n + e~ 


3.7% 


1.5% 


0.3% 


4% 


2% 


6.0% 


• K + e + e~ 


3.7% 


1.1% 


0.3% 


4% 


2% 


5.9% 




3.7% 


2.3% 


0.3% 


5% 


2% 


6.9% 


K + e+n~ 


3.7% 


1.6% 


0.3% 


4% 


2% 


6.0% 


K+n + e~ 


3.7% 


1.6% 


0.3% 


5% 


2% 


6.7% 


pe + e~ 


2.5% 


1.0% 


0.3% 


2% 


2% 


3.9% 




2.5% 


2.3% 


0.3% 


4% 


2% 


5.6% 




2.5% 


1.7% 


0.3% 


3% 


2% 


4.7% 


p/i + e~ 


2.5% 


1.7% 


0.3% 


3% 


2% 


4.7% 


■ ^ + <Pe+e- 


3.7% 


0.8% 


0.3% 


3% 


2% 


5.2% 




3.7% 


1.7% 


0.3% 


5% 


2% 


6.8% 


■ K + (f> e + e - 


3.7% 


0.8% 


0.3% 


3% 


2% 


5.2% 




3.7% 


1.6% 


0.3% 


5% 


2% 


6.7% 
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For the background shape assumption, the signal fits are repeated using a second order poly- 
nomial as the background PDF instead of the nominal first order polynomial. The result of the fit 
with the higher limit is used to quote the upper limit, unless variations of the signal PDF yield an 
even higher limit. 

In the normalization modes, the statistical uncertainties from the fits, the MC statistics and 
uncertainties from the signal and background shapes are all at or below 1%. The main uncertainty 
related to the normalization modes comes from how the efficiency is affected by sub-resonances 
in the decays. This is estimated to be about 3.5% for each of the modes and 1.1% for the 
A£ — ► pK~ir + mode. 

The requirement that the distance of closest approach between the two leptons be below 250 /im 
is not applied to the normalization mode and any MC-data difference can therefore not be expected 
to cancel in the branching ratio. Studies of the J/ip — > £ + £~ samples show efficiency differences 
between data and MC to be less than 0.3%. 

The efficiency of the particle identification has associated systematic uncertainties. We assign 
0.6% for each pion, 1.1% for each kaon, 1% for each electron and 2% for each muon. We do 
not apply a systematic uncertainty for the protons, since both the signal and the normalization 
mode contain a proton and the uncertainty therefore cancels. Uncertainties from the same types 
of particles are added linearly, while for different types they are added in quadrature. 

The efficiencies of the signal and normalization modes do not have the same dependence on p* . 
Varying the p* distribution used in the MC to better match the data changes the ratio of signal 
and normalization mode efficiencies by less than 2% for all decay modes. 

In the calculation of the signal efficiency we assume that the decays follow a three-body phase- 
space model. The selection efficiency has some dependence on where the decay lies in the Dalitz 
plane, so this assumption introduces a systematic uncertainty. Ignoring the regions we explicitly 
remove in the selection and the very high end of the m(£ + £'~) distribution, the efficiency typically 
varies less than 25% around the average as a function of m(£ + £'~). This model dependence is not 
included in the systematic uncertainty. 



5 PHYSICS RESULTS 

The invariant mass distributions of the normalization decay modes and the corresponding fits are 
shown in Fig. El The ir + K + K~ invariant mass distribution has both D + and Dj signals and is 
shown multiple times, because the D + and Df selections have different p* requirements and the 

— > ir + £ + £'~ selection requires a 7 from a D* + decay. Each selection requires its own fit in order 
to cancel most systematic uncertainties. The fitted signal yields are listed in Table 01 The table 
also lists the efficiencies estimated from signal MC. 

The invariant mass distributions for signal candidates in all 20 decay modes are shown in 
Figs.0HHl The yields obtained from unbinned likelihood fits are listed in Table|I]with statistical and 
systematic uncertainties. Only systematic uncertainties associated with the signal and background 
PDFs are included in the systematic uncertainty for the yields. The curves representing the fits 
are overlaid in the figures. No significant signals are seen and we calculate upper limits on the 

r(D+ .->%+e+i'-) v(d+ ^k+i+i'-) r(A + -*pe+e'-) <» 
branching ratios — „,v:_i. 7— t — , — :rr^ — and -^ttt — ^ — zr at 90% CL. For comparison 

r(D+ ) ^7r+0) r(D+ ) ^7r+0) r(A+^ P K-TT+) 

with previous measurements, the upper limits on the total branching fraction (BF) at 90% CL, 
calculated using Tabled are also given. 

The background from misidentified Df — > tt + tt +/ k~ decays in the Df — > 7r + /i + /i~ decay channel 
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Table 3: MC efficiency and fitted yields in data for the normalization modes. The ir + K + K~ 
distribution is fitted after the different selection criteria matching those for the corresponding 
signal modes are applied. 



Decay 


mode 


Select 


ion Mode 


N sig 


Efficiency 


D+ -> 


^ + <t>KK 




■ TT + /K+e+i'- 


25825±200 


(4.77±0.04)% 


Df-^ 


^ + <t>KK 


Df- 


> TT + £ + l'- 


23372±162 


(1.26±0.02)% 


Df^ 


^ + <t>KK 


Df- 


> K+t+a- 


71501±307 


(3.77±0.04)% 


At^ 


pK~7T + 




p£+£'- 


212664±1028 


(7.38±0.07)% 


D+ -> 




D+ - 


' K + <t>l+£- 


32920±329 


(6.12±0.05)% 




TT + (f)KK 


Df- 




115480±413 


(6.49±0.05)% 



is found to be 5 ± 2 times larger than expected from the weighted MC, while in all other muon 
modes no significant component of this background is visible. Cross checks of the MC did not 
lead to an explanation for the high level of this background in the — » 7r + /i + /x~ decay channel. 
Fitting the 7r + /i + /x~ candidate mass distribution with the misidentification component fixed to 
gives a signal yield of 0.2 events with little change in the upper limit. 

Figure |H1 summarizes the previously published and the new limits. For most decay modes, 
this analysis gives a significant improvement over the existing measurements. A recent CLEO 
measurement [7] is slightly more sensitive to D + — » ■n + e + e~ decays, and our sensitivity to D + — > 
7r + /x + /i~ and D + — > K + fj, + fi~ decays is worse than fixed target experiments [9] due to the low 
muon efficiency and high backgrounds. The D0 collaboration has recently presented a preliminary 
result [8] that improves the limit on D + — > Tr + /j, + (i~ decays to 4.7 x 10 -6 . 

The Dfs — > 7r + e + e - and — » 7r + jU + / _ t - candidates are shown in Fig. EH Signals are seen for 
all decays except for D + — * <t> u+ u~ ■ Table|5]gives the fit yields and eft — > £ + £" branching fractions 
calculated by normalizing to the Dfs — > tt^^kk decay modes. The significance is calculated from 
the change in likelihood with and without any signal. The observed branching fractions are in 
agreement with the world averages [11] of (2.98±0.04) x 10~ 4 for — * e + e~ and (2.85±0.19) x 10~ 4 
for cj) 
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Figure 3: Invariant mass distribution for ir + (f) K + K - candidates with a) D + — ► tt + £ + £' and 
D + -» K + £ + £'~ selection criteria, b) D+ -> tt + £ + £'~ selection criteria and c) D+ -> K + £ + £'- 
selection criteria, d) Invariant mass distribution for A£ — > pK~ir + candidates with the 7l+ — > 
p£ + £'~ selection criteria. The solid lines are the result of a fit to double- Gaussian signals and a 
second-order polynomial for the background. 
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Tt*e + |T invariant mass (GeV/c 2 ) invariant mass (GeV/c 2 ) 

Figure 4: Invariant mass distribution for D + — * ir + l + l'~ candidates. The solid lines are the results 
of the fits. The misidentified background component in the dimuon mode is shown as a dashed 
curve. 
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Figure 5: Invariant mass distribution for — > 7r + ^ + ^ /_ candidates. The solid lines are the results 
of the fits. The misidentified background component in the dimuon mode is shown as a dashed 
curve. 
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Figure 6: Invariant mass distribution for D + — ► ET + £ + £ /_ candidates. The solid lines are the 
results of the fits. The misidentified background component in the dimuon mode is shown as a 
dashed curve. 
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Figure 7: Invariant mass distribution for Df — > i^T + £ + £ /_ candidates. The solid lines are the 
results of the fits. The misidentified background component in the dimuon mode is shown as a 
dashed curve. 
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Figure 8: Invariant mass distribution for A£ p£ + £'~ candidates. The solid lines are the results 
of the fits. The misidentified background component in the dimuon mode is shown as a dashed 
curve. 
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Figure 9: Comparison of the branching fraction limits measured in this note with previously published measurements [7,9, 10]. 
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Figure 10: Invariant mass distribution for (left) D 
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and (right) D 
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candidates. The solid lines are the results of the fits. The misidentified background component in 
the dimuon mode is shown as a dashed curve. 



6 SUMMARY 



has been 



A search for the decay modes Df, -> -k + 1 + £'~ , D~£\ ->• K + £ + £'~ and A+ -> -pt 
performed using 288fb -1 of e+e - data. No signals are observed and we obtain upper limits on the 



branching ratios — 



-i+t'-) r(vf s) -*K+e+e'-) 



and 



r(A+^ P £+e-) 



between 10" 4 and 40 x 10" 



at 90% CL. This corresponds to limits on the branching fractions between 4 x 10 -6 and 4 x 10 -5 . 
These limits are calculated under the assumption of three-body phase-space decays; the efficiency 
varies by up to 25% as a function of dilepton invariant mass. For 17 out the 20 decay modes, the 
limits are an improvement over the existing measurements. 
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Table 4: Yields from fits to the candidates in the 20 Y+ — ► h + £ + £'~ decay modes. The first error is 
statistical and the second the systematic error on the yield. The third column shows the estimated 
signal efficiency. The fourth column shows the 90% CL upper limits on the branching ratios of 
the signal mode to the normalization mode. The last column shows the limits on the branching 
fraction for the signal modes at 90% CL. The upper limits include all systematic uncertainties. 

Yield BR (lO^ 4 ) BF (10~ 6 ) 

Decay mode (events) Efficiency (90% CL) (90% CL) 

D + ~ n + e + e - 24.0^-0+3.4 

d+ -> vr w- i-s±?g:i±i:S 

D+ -> TT+e+M" 4.1±i5:|±|:i 

n+ + -12 i + 15 - 5 + 3 - 2 

u — > 7T ^ e iz.i_ 14 8 _ 

d+ - vr w -9-4±5:2±?:2 

n+ _> 7T+P+//- 48 +4.7+o.8 

— > 7T e ^ 4 -°-3.9-0.3 

£>+ -> vr+^+e- 0.5j£°tJ-° 

D+ ^+ e +e- 5.9±f"l±g-| 

-> *W 2.9t?;°tg;? 

- i^ + e V- -3.41^+1-0 

D+ _> ^+ e + /U - -3.7±ti±i1 
D+ - -6.5±t:§±? : ? 

4- -> pe +e- o«»:: :i ;; : ) 

4+ - wV" 6.9l|?l°;i 
A+ -> pe +u- 0.2+2.9+0.5 
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< 24.4 


2.27% 


< 17.1 


< 10.8 


2.29% 


< 9.3 


< 5.9 


1 14% 


< 2 1 


< 7 6 


0.31% 


< 5.1 


< 18.5 


0.66% 


< 6.2 


< 22.3 


0.65% 


< 3.8 


< 13.9 


o.zlTo 


^00 


< 5.2 


0.75% 


< 22.2 


< 14.0 


1.64% 


< 5.7 


< 3.6 


1.64% 


< 5.9 


< 3.7 


2.81% 


< 1.8 


< 6.6 


0.68% 


< 7.1 


< 25.4 


1.40% 


< 1.5 


< 5.6 


1.40% 


< 1.0 


< 3.6 


4.11% 


< 0.7 


< 3.6 


0.67% 


< 8.1 


< 40.4 


1.19% 


< 1.8 


< 8.9 


1.18% 


< 1.5 


< 7.5 



Table 5: Yields from fits to the Dt^ — > tt + (/)£+£- candidates. The first error is statistical and 
the second is the systematic error on the yield. The third column is the branching fraction for 
4> — > £ + £~ calculated by normalizing to the D^ — > 7t + 4>kk decay modes. The fourth column is 
the significance of the signal. 



Decay mode 




Yield (events) 


Efficiency 


BF(cf> -► t^r) (10- 4 ) 


Significance 


D + -> 7T+0 e 4 




in n+7.1+0.4 
ly - u -5.9-2.9 


4.11% 


O Q + l.5+0.2 

°- y -1.2-0.6 


2.7<7 


D+ -> vr+c/y 




1 O+4.5+0.1 
i -°-3.7-0.7 


0.67% 


< 8.2 at 90% CL 


0.3<r 


D+ -> vr+0 e 4 


e~ 


en c+9.0+1.6 
ou -°-8.1-3.9 


1.19% 


9 0+0.5+0.2 

z -°-0.5-0.3 


6.8(j 






1 c q+6.1+1.3 
10-O_ 4 g_Q 3 


1.18% 


q c+1.4+0.4 
°-°-l.l-0.2 


3.2(7 
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